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We use shear reversal simulations to explore the rheology of dense, non-Brownian suspensions, 
resolving lubrication forces between neighbouring particles and modelling particle surface contacts. 
The transient stress response to an abrupt reversal of the direction of shear shows rate-independent, 
nonmonotonic behaviour, capturing the salient features of the corresponding classical experiments. 
Based on analyses of the hydrodynamic and particle contact stresses and related contact networks, 
we demonstrate distinct responses at small and large strains, associated with contact breakage 
and structural re-orientation, respectively, emphasising the importance of particle contacts. Conse¬ 
quently, the hydrodynamic and contact stresses evolve over disparate strain scales and with opposite 
trends, resulting in nonmonotonic behaviour when combined. We further elucidate the roles of par¬ 
ticle roughness and repulsion in determining the microstructure and hence the stress response at 
each scale. 


I. INTRODUCTION 


II. SIMULATION MODEL 


The flow behaviour of dense suspensions is strongly 
affected by details of the microstructure and interparti¬ 
cle forces [T]. Recent theoretical [5], experimental Hi] 
and computational mm work suggests that particle sur¬ 
face contacts make a major contribution to suspension 
rheology, though their precise role, and importance rel¬ 
ative to hydrodynamic interactions, is still debated. To 
this end, shear reversal experiments, in which the flow 
direction is suddenly reversed, prove to be an elegant 
means of estimating the contact stress (CS) contribution, 
while probing microstructural anisotropy. For example, 
it was shown that following a flow cessation period, the 
shear [7] (or normal 0) stress in a suspension of ~ 40 ^m 
polystyrene spheres reaches an “immediate” peak upon 
reversal, then evolves nonmonotonically over a strain of 
« 3 to its steady state. The large-strain evolution was 
attributed to microstructural realignment [7]; the initial 
stress peak was hypothesized to represent the hydrody¬ 
namic stress (HS), leading to a suggestion of the larger 
role played by particle contacts in denser suspensions [5| . 
It remains difficult to isolate the evolution of the CS and 
HS contributions and link the microstructural effect to 
the puzzling nonmonotonic behaviour. In the present pa¬ 
per, we reveal responses at two disparate strain scales, an 
elegent manifestation of the microfragile vs. macrofrag- 
ile distinction proposed by Cates et al. H- The small 
strain stress peak is shown to be a hydrodynamic re¬ 
sponse to surface-contact breakage, but is distinct from 
the steady state HS. The large strain scale is determined 
by microstructural reorientation, as predicted [Tj. The 
nonmonotonic behaviour is a combined effect of the CS 
and HS evolutions. We show that different surface char¬ 
acteristics control the stress response at different strain 
scales, meaning our two-scale explanation, and hence the 
micro- versus macro- fragility paradigm, can be applied 
usefully to a wide range of suspended systems. 


We solve the equations of motion numerically m 
for neutrally buoyant suspended non-Brownian particles, 
subject to forces and torques arising due to hydrodynam¬ 
ics and particle surface contact mi. For dense suspen¬ 
sions, in which the average surface separation between 
neighbouring particles becomes very small, the full hy¬ 
drodynamic resistance matrix [12] can be suitably ap¬ 
proximated by resolving pairwise, frame-invariant lubri¬ 
cation forces [13] . which diverge at contact and signifi¬ 
cantly exceed the long-range force components. Such a 
simplihcation has been proven to be effective in captur¬ 
ing the behaviour of dense suspensions Hmug. For an 
interaction between particles i and j, (with particle and 
fluid density p) the force and torque on particle i due to 
hydrodynamic lubrication can be expressed as 
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for particle diameter di, fluid viscosity rjf, particle trans¬ 
lational and rotational velocity vectors and cOi respec¬ 
tively, centre-to-centre unit vector n^- pointing from par¬ 
ticle j to i and identity tensor I, with the squeeze Qsq, 
shear Ugh and pump Upu resistance terms as derived by 
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The separation between particles i and j is calculated 
according to h = |rij| — for centre-to-centre vector 

r. We calculate the lubrication force when the interpar¬ 
ticle gap h is smaller than dmax = 0.05d (where d is 
the harmonic average particle diameter). An increasing 
body of evidence El 111 shows that direct particle-particle 
surface contacts can play a major role in suspension rhe¬ 
ology; indeed, simulations that strictly resolve lubrica¬ 
tion forces (treating particles as ideally hard and ideally 
smooth) [T7] have proven to be inadequate for captur¬ 
ing dense suspension rheology for cases where particle- 
particle contacts are presumed to be important. There¬ 
fore we truncate the lubrication divergence and regularize 
the contact singularity at a typical asperity length scale 
drnin (= O.OOld unless specified otherwise), i.e., setting 
h = dniin in the force calculation, when h < dmin- We 
use a value of r]f = 0.1 [viscosity unit: pd'^/t]. 

Mechanical contact occurs at /i < 0, giving normal re¬ 
pulsive and tangential forces described by a linear spring 
model and related through a Coulomb friction coefficient 
pLp (= 0.2 unless specified otherwise) [TT]. A linear (as 
opposed to Hertzian) spring is chosen for convenience, 
though we expect Hertzian results to lead to identical 
conclusions reagarding the respective roles of contacts 
and lubrication. The normal (F'^’") and tangential (F'^’*) 
contact force and torque F'^ are given by 
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The bulk stress tensor is calculated from the parti¬ 
cle force and velocity data. It is decomposed into con¬ 
tributions due to the hydrodynamic interaction and the 
particle-particle interaction, given by Eqs. [4a1and[4bl re¬ 
spectively, 

i i=jLj 

i i=jLj 

In the following discussion, we consider the shear com¬ 
ponents of the above stress tensors corresponding to the 
direction of the applied deformation, cr* and cr'^, as well as 
the mean of the diagonal components, namely the “pres¬ 
sures” and P'^. The hydrodynamic stress cr* is further 
decomposed in two ways. In the first, we isolate the con¬ 
tributions from normal forces (the squeeze Ugq terms) and 
tangential forces (the shear Ush and pump ap„ terms) [T3] 
as alarmed and respectively. In the second, we 

isolate contributions from opening and closing particle 
pairs (pairs for which dh/dt > 0 and dh/dt < 0 respec¬ 
tively), presented as cr^^pening and It is noted that 

^opening + fTclosing = and CT^rmal + f^Lngential = As¬ 
semblies of 5000 spheres are sufficiently large to achieve 
system size independence, and bidispersity with diame¬ 
ter ratio 1 : 1.4 prevents crystallisation m- Simulation 
results are ensemble averaged over 20 realizations with 
different initial particle configurations. We note that al¬ 
though the overlap is exceedingly small, typically of or¬ 
der 10 “^d in the Stokesian regime m, it can lead to 
qualitatively different rheology from that produced using 
the “ideal” hard-sphere model, as demonstrated later. 
The present technique produces results at the dense limit 
(solid volume fraction </> > 0.45) closely approximating 
those that would be obtained by fully resolving the hy¬ 
drodynamics (e.g. [ 20 ])) but assuming particles co-move 
with fluid at the mean flow level m, valid for shear flows. 
We verify this by incorporating an additional drag force, 
similar to [14] , which leads to a negligible increase in the 
calculated suspension viscosity. The particle assemblies 
are subjected to rate ( 7 ) controlled simple shear flow in 
a 3-dimensional periodic domain at constant (j) (= 0.54) 
and Stokes numbers St (= p^d^/pf) < inhibit¬ 

ing particle inertia. The suspension is first sheared from 
qt = —8 —>■ —2, reaching steady flow. No shear is ap¬ 
plied for qt = —2 —>■ 0. From yt = 0, the suspension is 
sheared in the opposite direction until a new steady state 
is obtained. 


(4a) 

(4b) 


III. STRESS AND MICROSTRUCTURE 
EVOLUTION 


for a collision between particles i and j with normal and 
tangential spring stiffnesses kn and kt respectively [A:„ = 
20000, unit: pd^/t^ and kt = (2/7)fc„], particle overlap 5 
and tangential displacement uy. 


The total stress (cr = cr“ -I- a’’) evolution. Fig l[a), 
is strikingly reminiscent of classical experiments [T] jS] 
M- Rate-independence is demonstrated by collapsing 
stress components with the respective steady state total 
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FIG. 1. (a) Evolution of stresses following reversal at jt = 0, 
showing total (cr, P), contact {a‘^, P‘^) and fluid (u*, P*) 
contributions, each scaled by the steady state total stress 
(ft, P). Multiple lines for cr illustrate j independence. In¬ 
set: Same data with logarithmic x-axis; (b) cr* components 
arising from (i) normal and tangential forces; (ii) opening 
and closing interparticle gaps. Inset: Evolution of the mean 
fluid film thickness {h) and scaled mean neighbouring-particle 
normal (I4i) / (Ki) and tangential (Vi) / (Ki) velocity magni¬ 
tudes. Embedded: particle-pair conhgurations corresponding 
to different times, (c) Evolution of coordination number Zc 
and shear fabric component Ai 2 - Inset: Coordination Zc and 
surface coordination Zh^^^ evolution, with logarithmic a;-axis. 


stress (T, for multiple j. Fig. [^c) shows the microstruc- 
tural evolution, characterised by a mechanical coordina¬ 
tion number Z^, the mean number of per particle contacts 
which support a contact stress greater than 10“® of the 
mean steady-state stress P, a surface coordination num¬ 
ber Zh„i„ counting all pairs with h < hmin and a fabric 
tensor HI], A = 2/(Zh_ A^) - il. 

Under shear flow, particles preferentially align along the 
compressive axis at 45° with the shear component of A, 
IA 12 I = 0.5 representing perfect alignment of all contacts 
and Ai 2 = 0 representing perfect isotropy. The evolu¬ 
tion of additional microstructure variables is shown in 
Appendix [Aj Fig|^ 


A. Steady flow and cessation 

In the steady state (jt < —2, Fig[^a)), the contact 
stress contribution is surprisingly large given the small 
contact overlaps, representing 60% of the total shear 
stress. The relative contribution is ^-dependent, e.g., 
at (/) = 0.47, we found cr° « 0.3cr [19]. We find Zc « 1.5 
and the shear component of A, A 12 « —0.01, indicating 
persistent mechanical contacts and an anisotropic net¬ 
work of lubrication films. In this condition, relative par¬ 
ticle motions are as illustrated in particle-pair diagram 
A, Fig[^b), a configuration that results in the mean rela¬ 
tive normal velocity of neighbouring particles (U„) being 
smaller than the mean relative tangential velocity (Ut), 
highlighted in Fig [^b) Inset, which gives these quanti¬ 
ties scaled by the steady state value of (U„), (Ki). This 
leads to comparable normal and tangential lubrication 
forces (and corresponding stress contributions cr^^ormai 
and as decomposed in Fig. j^b)), in spite of 

the order of magnitude difference expected from their 
respective 1/h versus \n(l/h) dependence. The average 
number of particle pairs moving together or apart is equal 
at steady state, as required to satisfy the constant vol¬ 
ume constraint, resulting in constant mean lubrication 
film thickness (h), Fig. Sb) Ins^ equal stress contribu¬ 
tions cr(,pe„i_ and CTeiosing, Fig. mb), and P dominated 
by Pj Fig[5a). 

Upon flow cessation (jt = —2), the contact stress re¬ 
laxes together with the hydrodynamics stress, suggest¬ 
ing that caution should be exercised when interpreting 
the “instantaneous” stress loss as entirely hydrodynamic 
in such experiments [T| [23j . Correspondingly, Zc drops 
to zero in the relaxation period, though a small portion 
of weak contacts relax more slowly due to confinement 
and fluid overdamping. The shear-induced anisotropic 
microstructure pertaining to hydrodynamics, however, 
remains intact throughout the relaxation period, evi¬ 
denced by constant A 12 and (Fig [^c) and Inset, 

^ 0) = Zh^.^(jt = 4)), implying the steady 
state HS can be recovered instantaneously (with oppo¬ 
site sign) upon shear reversal. 
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B. Shear reversal: micro- and macro-strain 
responses 

Indeed, a’’ does resume it’s steady state magnitude 
upon reversal for strains < 10“'*, Fig.ma) Inset. It then 
surges to a significant peak, around 50% greater than the 
steady value, at strain 10“^, sustaining until about 10“^ 
where it starts to subside. Resumption of the steady 
value followed by a demonstrable peak is also observed 
for P* over the same strain scale. We attribute this small- 
strain surge, the manifestation of a microfragile response 
[5], to the pulling apart of particle surfaces at the /imin 
(= 10“^(i) scale due to the new (reversed) load being in¬ 
compatible with the present microstructural alignment. 
This is clearly demonstrated by the coincident decrease of 
which reaches a minimum near 10“^. The mech¬ 
anism is further evidenced by the significantly greater 
(Vn) than under steady flow, the dominance of 
and cr(,pg„ijjg, and the tensile nature of This is further 
illustrated by the evolution of the h distribution, given in 
Appendix Relative particle motions during this time 
are illustrated in particle-pair diagram B, Fig.j^b). Such 
microfragile events in, e.g., a dry granular system, would 
be subtle to detect or difficult to distinguish from the 
macroscopic process. These events in dense suspensions, 
nonhydrodynamic in nature, however, lead to the spec¬ 
tacular hydodynamic responses, which have been mea¬ 
sured robustly in experiments mm- We note that a 
microfragile hydrodynamic response is absent in Stoke¬ 
sian Dynamics simulations of shear reversal that strictly 
inhibit fluid films smaller than O.Old [23], strengthening 
the argument for direct surface contacts in addition to 
hydrodynamics, as a crucial contributor to the rheology 

observed by mm- 

The subsequent building up of after a strain of 

0.01 is coupled to re-orientation of the microstructural 
anisotropy A 12 , corresponding to macrofragile evolution 
at a larger strain scale of order unity. The initial subsi¬ 
dence of a’’ from its peak until jt « 0.5 (while A 12 < 0), 
corresponds to a net opening of lubrication films (see {h) 
and (Topening, Fig[^b)), consistent with the leading 1/h 
dependence of the lubrication forces, combined with a re¬ 
duction in (Vn)- At larger strains, a new contact network 
establishes in the now-compatible, oppositely aligned, 
compressive direction (evidenced by A 12 > 0) with net 
repulsive lubrication forces during 0.5 < jt <2, restoring 
(h) to its steady value thereby producing positive and 
a marginally dominant The consequent mean 

relative particle motion is highlighted in paritcle-pair di¬ 
agram C, Fig Sb). Although (T* evolves continuously 
during this large-scale period, the responsible mechanism 
therefore switches as the anisotropy changes sign. The 
stress presented by |23| has a comparable macrofragile 
evolution, but exhibits nonmonotonic behaviour due to 
the absence of a microfragile response. 

The contact stresses (ct'’, P'’) follow a similar 

macrofragile evolution, their associated microfragile con¬ 
tact breakage having occurred at flow cessation as dis¬ 


cussed. The stress evolution is closely correlated with 
the building up of the mechanical coordination number 
Z^, which occurs on a similar strain scale as the fabric 
reorientation described above, as illustrated in Fig[^c). 
The separation of scales in the evolution of the contact 
stress and the hydrodynamic peak ensures dominance 
of the hydrodynamic stress at small strains after rever¬ 
sal, an assertion made in BIH], though overlooking the 
microfragile hydrodynamic response. Combining the in¬ 
creasing cr^ with the decreasing cr* at yt > 0.01 gives rise 
to the nonmonotonic total stress, meaning microfragility 
in the hydrodynamic response is crucial in capturing the 
experimental behaviour. 

The above analysis sheds light on the two-scale nature 
of the stress evolution, linked to configurational change 
at small strains and anisotropy re-orientation at large 
strains. The importance of particle contacts in achieving 
the nonmonotonic stress response naturally leads to the 
question of the sensitivity of the evolution at each scale 
to particle interactions and surface properties. 


IV. ROLE OF PARTICLE PROPERTIES 

In order to test the applicability of the above described 
mechanism to a wide range of particle systems, we ad¬ 
dress two factors pertaining to well studied suspensions, 
namely surface roughness and stabilising repulsion. For 
suspensions of large particles (e.g., d > 10 /rm), such as 
the 40 — 50^m polystyrene spheres suspended in density 
matched silicon oils studied in |7], surface roughness is 
perhaps the more relevant factor; for those of small par¬ 
ticles (e.g., d < 10 fj.m) steric or electrostatic repulsion 
may give well-defined repulsive forces. 


A. Surface roughness 

Surface roughness is represented numerically by an as¬ 
perity length scale (by means of dmin) and the friction 
coefficient yip. dmin contributes to the strain scale of the 
microfragile HS response and should also affect the HS 
magnitude. We explore such effects by varying dmin be¬ 
tween 10“*d and 10“^d, considering the physical size of 
surface asperities and bounded numerically by the singu¬ 
larity and the overdamping requirement at the lower and 
upper limits, respectively. The resulting cr* scaled by the 
steady state total stress d, following a reversal at yt = 0, 
is plotted against strain on a log-linear scale in Fig.j^a). 
The strain scale of the microfragile peak decreases rather 
linearly with decreasing dmin in the O.Old to O.OOld range, 
but saturates approaching 10“*d. We verified that the 
saturation is not due to inertial effects, but is perhaps 
due to the nonlinear coupling between particle configu¬ 
ration and dynamics. Decreasing dmin also significantly 
increases the peak magnitude, but only weakly affects 
the macrofragile HS response, the evolution of which in¬ 
terestingly collapses relative to d. The surface roughness 
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FIG. 2. Stress response to reversal as a function of (a) asper¬ 
ity length scale /imin, with the numerical values of /imin given 
in the legend in units of d; (b) friction coefficient fip, with 
the numerical values of fip given; (c) repulsive force magni¬ 
tude |F''|, quantified as the relative magnitude of the repulsive 
stress d’' to the contact stress d-"^. 


effect of /i„iin thus controls the microfragile, but not the 
macrofragile HS response. 

In a manner often employed in dry granular studies 
(see e.g. [H]), we vary the friction coefficient /ip incre¬ 
mentally between 0 and 1, exploring particle surfaces 
from ideally frictionless to very frictional. The total shear 
stress evolution is given on a linear cc-scale in Fig. [^b). 
On the contrary to the effects of varying hmin, the mi¬ 
crofragile response is largely insensitive to fXp, which is 
unsurprising given the small scale stress is dominated by 
attractive lubrication forces. The invariance of the hydro- 
dynamic stress gives strong support to the central role of 


particle contacts in achieving the very different viscosities 
observed in such systems. The stresses differ hugely at 
larger strains, however, indicating that the /J,p-dependent 
contact stress is important during 0.3 < yt < 2, coincid¬ 
ing with recovery of Zc to its steady value. The increase 
of the contact stress with increasing friction can be un¬ 
derstood from the increase of tangential contact forces 
and the decreased departure from the jamming volume 
fraction (pc |2fij , which is known to decrease as friction in¬ 
creases [221 Ei]. The latter effect is also consistent with 
the experimental observation that the peak immediately 
after reversal becomes lower relative to the steady state 
stress when increasing volume fraction |^. The inter¬ 
particle friction thus mainly affects the large scale mi¬ 
crostructure and contact stress and hence the macrofrag¬ 
ile response. In reality, /ip and /imin are probably si¬ 
multaneously coupled to the surface roughness variation, 
though the combined effect may be deduced from the 
present separate analyses, exploiting the marked separa¬ 
tion of scales associated with our two-scale description. 


B. Surface stabilisation 

We next probe the effect of a generic stabilising re¬ 
pulsive force, extending the above analysis to consider 
particles in the size range d < lO/rm. It is assumed, 
based on previous simulation results 0161, that a static, 
short range, normal repulsive potential is sufficient to 
capture the essence of a stabilising mechanism such as 
electrostatic repulsion or a grafted polymer hair coat¬ 
ing. Enhanced dissipation in the lubrication forces, a 
phenomena described by [55], is neglected for simplicity. 
A generic form of the repulsive force model derived by 
Fredrickson et al. [21] is used, F’' = A: , where k 

is some constant that encapsulates (among other things) 
the chemical properties of the hairs and their density on 
the surface, essentially quantifying the “strength” of the 
static repulsion. We apply the same singularity regular- 
isation as in the lubrication model, and the same values 
^min, hjnax- Coupling to the mechanical contact model 
is as before. The total shear stress response to reversal is 
given in Fig. j^c), for k spanning 2 orders of magnitude 
(quantified by the relative magnitude of the steady-state 
repulsive ct” and contact (t° stresses). As expected for 
small fc, the additional static repulsion is insufficient to 
separate particles, so the stress response closely resem¬ 
bles that for the base case in Fig. [^a). As k (or a'" ja'^) 
is increased, we note that while the large strain scale for 
the evolution appears to be unchanged, the steady value 
of (T decreases. This is attributed to increasing inhibition 
to mechanical contacts (for which ft, < 0) as the repul¬ 
sion becomes stronger. We note that this trend is valid 
when (t“ is of comparable magnitude to o'”. For very 
large an opposite trend is observed 0 due to a 

shear thinning mechanism—the polymer hair length can 
begin to contribute to an effectively larger total particle 
diameter, leading to a higher effective volume fraction 















6 



) -- 1 - 

-6 -4 -2 0 2 4 


it 

(a) 



it 

(b) 


FIG. 3. (a) Response to shear reversal, showing the total 

a, contact cr^, fluid cr* and repulsive stress a'~ contributions, 
each scaled by the steady state total stress tj; Right Inset: 
Same data with logarithmic a:-axis. Middle Inset: Evolu¬ 
tion of mean fluid film thickness {h) and surface coordination 
Zhmin ■ (b) Analogous result for very short relaxation period. 


and therefore a higher shear stress, as explained in detail 
by [6] and references therein. This then leads to shear 
thinning behaviour with reducing k, rather than with in¬ 
creasing k as we observe here. For small strains after 
reversal, we observe a marked loss of the microfragile 
stress peak as k is increased. 

To further understand this loss, we present the full 
evolution of shear stress contributions for large k and a 
flow cessation period sufficient to relax to steady state. 
Fig. I^a), with (h) increasing from around O.Ofld to O.OSd 
and decreasing modestly (middle Inset). To further 
characterise the relaxation period, we provide the associ¬ 
ated h distributions in Appendix [Bj under the action of 
the repulsive potential. Upon reversal, some remaining 
timin contacts are opened, resulting in a decrease 

and a HS response over a 10“^ strain scale (middle and 
right Insets, respectively), consistent with the microfrag¬ 
ile response in Fig. In this sense, a microfragile HS 
response still occurs; though it starts from a “loosened” 


microstructure, producing a HS lower than its steady 
value, rather than the surge noted previously. Follow¬ 
ing this reasoning, a HS peak would be recovered if the 
relaxation period were shortened sufficiently to disallow 
any increase in (h). We verify this in Fig [^b) using a 
very short relaxation period it = —O.OI —0. A HS surge 
of about 100% of its steady value is observed, although it 
does not result in an appreciable peak in the total stress 
since the HS contribution is small. The repulsion and 
contact stresses similarly follow a macrofragile evolution. 
This again creates a strain window for the HS to be mea¬ 
sured separately from other components. In short, the 
repulsive force magnitude together with the associated 
scales provides extra control over particle configurations 
and hence the stress response. The two-scale evolution 
concept is, however, still robustly helpful in understand¬ 
ing this more complicated behaviour. 


V. CONCLUDING REMARKS 

We provide a robust characterisation of a two-scale 
response to shear reversal in dense suspensions, that is 
highly reminiscent of the micro- versus macro- fragility 
proposed by Cates et al. [5]. Links are established be¬ 
tween stress responses at small and large strains with mi¬ 
crofragile contact breakage and macrofragile microstruc- 
tural (re-)building respectively, resolving the hitherto un¬ 
explained nonmonotonic stress evolution following shear 
reversal. This substantiates the emerging understanding 
about the importance of particle contacts in suspension 
rheology - they not only provide a significant contact 
stress at steady state, but also give rise to a pronounced 
small strain transient hydrodynamic response. This un¬ 
derstanding provides a sound theoretical framework from 
which to formulate constitutive models with appropriate 
two-scale characteristics, and previous attempts at such 
models [30] might be revised to correctly link the stress 
and microstructure at each scale. The evidence that dif¬ 
ferent microstructural features control the contact and 
hydrodynamic stresses respectively and in an analogous 
way to that in dense granular flows |22|, supports fur¬ 
ther unification of dense suspension and granular rheol¬ 
ogy extending from steady m to unsteady state. The 
findings on surface features and interactions also open 
doors to either devising new experiments and protocols, 
e.g., varying relaxation time, to characterise particle sur¬ 
face properties and stress contributions; or designing new 
particles, e.g., with different grafted polymer hairs, to re¬ 
alise certain desired rheological properties. 


ACKNOWLEDGEMENTS 

This work is funded by the Engineering and Physical 
Sciences Research Council (EPSRC) UK and Johnson 
Matthey through a CASE studentship award. The au¬ 
thors would like to thank M. E. Cates and M. Hermes 































7 



7 

6.8 

I 

6.6 



(c) 


(d) 


FIG. 4. Further microstructural quantities to support Figure 
1(a). (a)-(b): Microstructural evolution evaluated at contact, 
to support Zc.\ (c)-(d) Microstructural evolution evaluated at 
lubrication dutoff, to support A\i. Definitions of Z and Acj ^2 
given in Appendix]^ Zc and A 12 as defined in main article. 


for their critical reading of this manuscript, and B. Guy, 
W. C. K. Boon, P. McGuire, M. Marigo, H. Stitt, H. Xu 
and J. Y. Ooi for helpful discussions. 


than 10 ® of the mean steady-state stress P, consistent 
with the definition of Zc given in the main article. 


Appendix B: Distributions of h 

We plot the distribution of the particle-particle sep¬ 
aration length h for the simulations in Fig. 1(a) and 
Fig. 3 of the main article. It is noted that in the case 
with significant polymer hair repulsion, there remains a 
peak in PDF{h) at very small h. This is consistent with 
the corresponding stress evolution, which demonstrates 
that there is still a non-negligible contribution from di¬ 
rect particle-particle contacts, cr'^. 



(a) (b) 


Appendix A: Further microstructural quantities for 
Figure 1 

Additional microscopic quantities to back up the find¬ 
ings reported in Figure 1(a) of the main article are pre¬ 
sented in Figure]^ of this document. Z is a lubrication 
contact number, counting all pairs with h < hmax- All 
such pairs contribute to the fabric component A 12 . 
is the shear component of the mechanical fabric tensor, 
which we define as Ac = 2f{ZcN) J2h<o ~ 

omitting those pairs which support a contact stress less 


FIG. 5. Distribution of h for the case in (a) Figure la and 
(b) Fignre 3 of the article, focussing on the range for which 
lubrication forces are calculated. 

In addition, we find that the peak in PDF{h) at small 
h remains even after the relaxation period. We attribute 
this somewhat counter-intuitive finding to the repulsive 
force magnitude and cut-off scales and confinement ef¬ 
fects. It is the subsequent opening of these remaining 
small h particle pairs that is responsible for the very rapid 
evolution of cr* reported in Fig. 3 of the main article. 
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